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Minijet production in ultra-relativistic heavy-ion collisions is discussed by taking semi-hard 
parton rescatterings explicitly into account. At both RHIC and LHC energies we find sizable 
effects on global characteristics of the nuclear collision like the initial multiplicity and trans- 
verse energy of the minijets. The dependence of these quantities on the cutoff that separates 
soft and hard interactions becomes much smoother after the inclusion of the rescatterings. 
This allows to define an energy and centrality dependent saturation cutoff and to push per- 
turbative computations to rather low values of the cutoff. As an application we compute the 
charged multiplicity at mid rapidity and compare it to the recent RHIC data. 

In nuclear collisions at ultra-relativistic energies the initial production mechanism is the 
liberation of a great number of partons from the nuclear wave-functions in a very short time 
from the beginning of the collision. Then, this dense system of partons, also called minijet 
plasma, will then evolve, possibly thermalizing and giving rise to the quark gluon plasma. The 
signals of its creation are very sensitive to the minijet plasma initial conditions, that in turn 
can be related to final state bulk observables like the charged particle multiplicity and their 
transverse energy. On the other hand, at RHIC and LHC an increasingly large fraction of the 
event is expected to be due to hard and semi-hard interactions, so that we can try to gain control 
on the initial conditions by using perturbative QCD. 

The usual way to do this is to take into account the Eikonalized minijet cross-section a m j = 
j d 2 (3 (1 — exp [— ct/Zab (/?)]), where f3 is the nuclear impact parameter, Tab(P) = / d 2 brA(b — 
(5)tb{o), oj = J dxdx' G(x)(Jh{xx')G(x'), G is the parton distribution function and t^b) are the 
nuclear thickness functions; <th{xx') is the perturbative parton-parton cross section defined with 
an infrared cutoff p C ut, that separates soft and semi-hard scatterings; here and in the following 
the flavour indices are suppressed for simplicity. In this way we are describing different parton 
pairs interacting in different points in the plane transverse to the beam; we will call these events 
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Figure 1: a) An example of disconnected collisions, b) Rescatterings. 

disconnected collisions, see Fig. |l]a. In each collision two back-to-back jets are produced, so that 
the minijet multiplicity is easily seen to be 



N^((3) = 2N col ((3)=2ajT AB (f3) 



dxdx' 'd 2 bT a(x , b — (3)o h{xx')T b(x' , b) 



(1) 



where Ta(x, b) = G(x)ta(P). 

There are two problems with this expression. First, it diverges when p cu t— ^0 due to the 
inverse power divergence of the perturbative cross-section an. Second, at very high energy and 
large atomic numbers a typical projectile parton is traversing a dense target nucleus, so that the 
chance of scattering more than once may become non negligible. To have a quantitative feeling 
on this, one can take the average number of parton-parton collisions in ([!]), drop the integrations 
over x and b, divide by the number of incoming partons, Ta, and obtain the average number 
of collisions per incoming parton, (N scat (x,b)} = J dx' 'cth(xx')T b(x' ,b). For example, at LHC 
with a cutoff p cu t = 2 GeV a parton in the centralrapidity region suffers on average two to three 
collisions over most of the target transverse area El. In conclusion, parton rescatterings must be 
taken into account in the description of the collision dynamics. 

The simplest rescattering event is a three-parton elastic scattering. It can be shown! that 
in the limit t/s— >0 and for infinite number of colours this cross-section factorizes in terms 
of two two-body elastic cross-sections. Therefore the process can be pictured as two succes- 
sive semi-hard scatterings against two different target partons (see Fig. ||b). We can try to 
extend this picture to an arbitrary number of scatterings in the following way i. We con- 
sider the multi-parton distributions D^(x,b), i.e. the probability densities of having n par- 
tons in the A nucleus with fractional momenta x% and transverse coordinates b{. If we ne- 
glect correlations inside the nuclei the multi-parton distributions are Poissonian: D^(x,b) = 
■^jTa(xi, b\) . . . TA(x n , b n ) exp[— / dxd 2 bTA(x, &)]. Then, we consider the probability that at 
least one semi-hard interaction occurs between n partons from A and m partons from B: 



Pnm = 1 — Iir=i ITj=i(l ~~ where &ij = aH(xiXj)5^ (6j — bj) is the probability of a semi-hard 



scattering between the i-th parton from A and the j-th from B. One can then obtain the mini- 
jet cross-section by convoluting the multi-parton distributions and the semi-hard probability: 
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Dg. In this way we are describing both the disconnected collisions and 



the rescatterings (see Fig. |l|b). The minijet multiplicity is no longer related in a simple way 
to the number of collisions, however if we choose a fixed scale Q = p cu t for the distribution 
functions and the running coupling constant, we can obtain a simple expression for itu: 



N jet ((3) = J dxd 2 bT A (x,b- 0) (l - e -Ux'a H (xx')r B (x',b)^ + { A ^ £j 
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We interpret it as the density of projectile partons Ta times the probability of having at least one 
semi-hard scattering, and the exponent in (^) as the opacity of the target. As p cu t goes to infinity 
the target opacity becomes smaller and smaller, therefore the projectile is traversing a very dilute 
target, the chance of scattering twice or more becomes negligible and the minijet multiplicity 
approaches the Eikonal estimate. On the contrary, as the cutoff decreases the opacity increases, 
the target becomes almost black and the multiplicity reaches a finite limit: roughly speaking 
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Figure 2: a) Minijet multiplicity dependence on the cutoff. Solid lines include rescatterings, dashed lines are 
the Eikonal estimate, b) Saturation cutoff dependence on the cm. energy compared to the EKRT model, 
c) Saturation happens when the solid line crosses the long-dashed one, see Eq.(^). 



This basically expresses the fact that a collision cannot free more partons than were present in 
the incoming nuclear wave-function. The effect of the rescatterings is shown in Fig.^a, where the 
minijet multiplicity tends to saturate below 2 GeV at LHC and 1 GeV at RHIC. Note that the 
dependence on the cutoff has become much smoother than in the Eikonal estimate, making the 
choice of the cutoff less critical. An analogous behaviour is observed also in the initial transverse 
energy and for the A;-factor dependence of these quantities §. 

It is possible to show that in (^) we are taking into account not only the jets that had n > 1 
semi-hard scatterings, but also those that had at least one semi-hard scattering but also any 
number of soft ones El. Indeed one semi-hard scattering is enough to shadow completely the soft 
ones, so that in (^) the soft cross-section does not appear. In this way, we are neglecting only 
the minijets that had just soft scatterings. This, together with (||), allows to define a saturation 
cutoff p S at'- we require that the number of jets be near the limit and we call p sa t the cutoff at 
which this happens; for example we can require 



Njet(Pcut = Psat) = 80% N Um (p cut = p sat ) . (4) 

Finally, we compute the initial conditions at this cutoff: iV^jf = Nj et (p cu t = Psat)- 

Now, we have to check if p sa t is far enough form Aqqd, in which case we can trust these 
perturbative computations and say that the semi-hard minijets saturate the initial production 
mechanism. At RHIC energies this method is at the border of applicability, and at LHC it 
should be reliably applicable (Fig.|2|b). The percentage chosen in (Q) is an arbitrary parameter, 
however it happens that p sa t sits where the plateau in the minijet multiplicity begins (Fig^p), 
so that choosing a different value does not change too much the results. 

With these tools we can study the charged particle multiplicity. We assume isentropic 
expansion of the minijet plasma in the central rapidity region and direct proportionality between 
minijets and final hadrons 0, N ch (^/s, (3) = 0.9 | N?^(^/s, j3) , and compute the number of 
participant nucleons at fixed energy and centrality in a Glauber model: N part = J d 2 bTA(b — 
(5) (1 — exp[— o"pp(- v /s)tb(6)]), where a pp is the inelastic proton-proton cross-section. To study 
non-central collisions with Wood-Saxon thickness functions, we have to enforce the collision 
geometry by assigning the nuclei a radius R and by computing p sat inside the overlap area 
defined by this radius. The reason is that in its periphery a nucleus is very dilute, therefore in 
that region it makes no sense to impose that it becomes black. The results are shown in Fig. 
|3]a as a function of the cm. energy and of the number of participants. 

For central collisions the dependence on the energy is strikingly similar to the EKRT sat- 
uration model 0, however in our case the saturation is determined dynamically and is related 
to the blackness of the target, moreover no shadowing corrections have been included in the 
parton distribution functions. For R equal to the Woods-Saxon radius, Rws, the slope of the 




VS (GeV) Npart 

Figure 3: Charged particle density at mid-rapidity per participant nucleon pair. 
Left: Center of mass energy dependence. Right: Centrality dependence. 



centrality dependence at yfs = 130 GeV compares very well with the experimental datau down 
to 180 participants; it underestimates phenix data but is not bad compared to brahms and 
phobos. Remember also that we are neglecting ah initio up to 20% of soft partons in the initial 
conditions computation, and that the saturation is not expected to work well for peripheral col- 
lisions. At this energy the dependence on the radius R is rather strong, but it decreases already 
at ^/s = 200 GeV, and almost disappears at LHC energies. Note also that the slope of the curve 
increases with the energy signaling an increasing role played by purely semi-hard events. This 
is in sharp contrast with the EKRT model, where an almost flat behaviour is found. 

In summary, the inclusion of the rescatterings in the dynamics of the collision allows to 
define global observables also at very low values of the cutoff, in a region of interface between 
perturbative and non-perturbative interactions. The minijet multiplicity tends to saturate and 
to reach a limit given by the number of partons in the initial nuclear wave functions. This allows 
to define a dynamical saturation cutoff at which the initial conditions can be estimated. The 
results compare satisfactorily to RHIC data at = 130 GeV, where the model is at the edge 
of its applicability. As the energy increases it should be more and more reliable. 
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